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COLLISIONLESS  COUPLING  IN  THE  AMPTE  ARTIFICIAL  COMET 

I .  INTRODUCTION 

The  in  situ  measurements  of  the  plasma  parameters  from  the  AMPTE 
(Active  Magnetospheric  Particle  Tracer  Explorers)  artificial  comet  releases 
(Valenzuela  et  al.,  1986;  Haerendel  et  al.,  1986)  in  the  solar  wind  provide 
us  with  a  unique  set  of  data  to  test  the  available  theories  on  the  subject 
of  collisionless  coupling  of  magnetized  plasma  streams  under  high  Mach 
number  conditions  (MA  >>  1)  (Haerendel  et  al.,  1986).  The  subject  is  very 
opportune  since  it  is  the  controlling  factor  that  determines  the  momentum 
coupling  process  occurring  in  the  interaction  of  the  solar  wind  with  the 
cometary  plasma  generated  by  the  ionization  of  the  neutral  coma.  The 
overall  comet  structure,  the  applicability  of  MHD  or  kinetic  models,  the 
presence  or  absence  of  a  cometary  shock  and  the  type  of  the  resulting 
ionopause  (Mendis  and  Houpis,  1982;  Ip  and  Axford,  1982;  Fedder  et  al., 
1986;  Sagdeev  et  al.,  1986)  depends  critically  on  the  wave-particle 
processes  producing  the  momentum  coupling,  the  thermal izat ion,  and  the 
isotropization  in  the  interaction.  It  is  the  purpose  of  the  present  note 
to  compare  the  AMPTE  in  situ  observations  of  the  plasma  parameters  and  wave 
signatures  (Gurnett  et  al.,  1985;  Haerendel  et  al.,  1986)  with  the 

theoretical  concepts  currently  applied  to  the  high  Mach  number  interaction 
problem.  In  the  next  section  we  present  a  brief  description  of  the 

experiment  as  well  as  the  relevant  data.  Section  III  reviews  the 

theoretical  models  for  coupling.  Section  IV  compares  the  AMPTE  data  with 

the  observations. 

II.  OBSERVATIONAL  RESULTS 

During  the  artificial  comet  experiment  on  Dec  27,  1984  two  canisters 
of  barium  were  released  from  the  IRM  (Ion  Release  Module)  spacecraft.  The 
release  was  on  the  morning  side  of  the  earth  at  a  geocentric  radial  dis¬ 
tance  of  roughly  17  earth  radii.  The  canisters  were  exploded  simultan¬ 
eously  at  a  distance  -  1  km  from  the  spacecraft  at  1232  UT.  The  explosion 
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produced  an  expanding  barium  cloud  which  was  rapidly  ionized  by  solar  UV 
( x ^  =  28  sec  where  is  the  pintoionization  time).  The  interaction  of  the 
solar  wind,  which  was  flowing  at  -  550  km/sec,  with  the  ionized  barium 
cloud  was  recorded  by  the  IRM  instruments  inside  the  cloud  and  the  magnetic 
cavity  that  was  created.  The  UKS  spacecraft,  located  -  170  km  away,  was 
outside  the  magnetic  cavity  and  measured  magnetic  disturbances  and  particle 
fluxes  generated  by  the  interaction.  A  schematic  of  the  spacecraft 
positions  and  the  magnetic  field  structure  is  shown  in  Fig.  1.  We  comment 
that  in  addition  to  the  in  situ  observations  made  by  the  IRM  and  UKS 
spacecraft,  ground  based  and  airborne  optical  data  were  also  obtained 
during  the  release.  Such  data  provides  important  information  on  the 
macroscopic  behavior  of  the  solar  wind  -  barium  interaction.  However,  for 
the  purposes  of  the  present  study  this  data  is  not  used  since  we  are 
concerned  with  the  detailed  evolution  of  plasma  and  field  quantities  which 
is  not  provided  by  optical  data. 

We  present  a  set  of  plasma  and  field  measurements  which  highlight  the 
dynamic  interaction  of  the  solar  wind  and  the  barium  cloud  in  Figs.  2-4. 
In  Fig.  2  we  show  measurements  of  the  electron  density  from  15  eV  to  30  eV 
[ne  (Fig.  2a)],  the  flow  velocity  of  the  solar  wind  protons  measured  in  the 
range  20  eV  to  40  keV  in  GSE  coordinates  [Vpx  (Fig.  2b)  and  Vpz  (Fig.  2c)], 
arid  the  magnitude  of  the  interplanetary  magnetic  field  [B  (Fig.  2d)] 
(Haerendel  et  al.,  1986).  We  note  that  the  magnetic  field  is  in  the  y- 
direction  and  that  the  electron  density  does  not  include  the  cold  electrons 
associated  with  the  barium  ions.  The  magnetic  cavity  generated  by  the 
barium  ions  is  clearly  seen  during  the  time  UT  12:32:02  and  UT  12:33:15. 
Figure  3  shows  the  low  frequency  electric  field  measurements  as  a  function 
of  time  and  frequency  (Gurnett  et  al.,  1985),  the  magnitude  of  the  magnetic 
field,  and  the  barium  ion  density.  Figure  4  displays  the  electric  field 
spectrum  upstream  of  the  ion  cloud  at  the  time  of  maximum  intensity  (UT 
12:34:27)  (Gurnett  et  al.,  1985).  The  electrostatic  waves  reached  ampli¬ 
tudes  in  excess  of  140  mV/m.  Since  the  electron  density  in  Fig.  2  does  not 
include  the  cold  elctrons,  for  times  prior  to  12:34:30  UT  we  will  use  the 
number  density  from  Fig.  3,  which  is  based  upon  wave  emissions  at  the 
plasma  frequency.  For  the  time  period  12:34:30  UT  -  12:35:30  UT  the  values 
of  the  electron,  proton,  and  barium  densities  are  not  well-known.  Finally, 
for  times  later  than  12:35:30  UT  we  use  the  data  presented  in  Fig.  2. 
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On  the  basis  of  these  measurements  the  following  picture  describing 

the  spatial  evolution  of  the  coupling  between  the  solar  wind  protons  and 

the  barium  emerges.  For  times  later  than  UT  12:36:20  the  solar  wind 

parameters  correspond  to  the  ambient  conditions  (n  =  2  cm  V  * 

e  px 

550  km/sec,  Tg  *  2  x  10  °K  and  B  »  10  Y ) .  The  interaction  between  the 
solar  wind  protons  and  the  barium  cloud  starts  at  the  point  marked  1 ,  which 
corresponds  to  UT  12:36:20,  as  seen  from  the  initiation  of  the  slowing  down 
of  the  solar  wind  (Fig.  2b).  During  the  time  period  UT  12:36:20  -  UT 
12:35:15  the  initial  slowing  down  rate  is  relatively  weak  and  is  accom¬ 
panied  by  moderate  wave  activity  (Fig.  3),  density  and  magnetic  field 
compression  (Fig.  2a, d).  There  is  no  plasma  flow  in  the  z-direction  (Fig. 
2c).  This  time  coincides  with  the  time  that  the  UK  spacecraft  records 
flux ss  of  hot  electrons  (>  100  eV).  A  much  stronger  slowing  down  rate  is 
observed  between  UT  12:35:15  and  UT  12:34:25,  accompanied  by  strong 
electrostatic  wave  activity  near  the  local  proton  lower  hybrid  frequency. 
For  B  *  30  -  90Y,  n  =  2  -  10  cm  and  na  3  2  -  120  cm-^  we  note  that  the 

p 

proton  lower  hybrid  frequency  is  fHp  3  10  -  30  Hz.  The  magnetic  field, 
density,  and  the  temperature  continue  to  increase,  while  the  value  of  Vp2 
remains  relatively  unchanged.  This  continues  until  UT  12:34:25  which  is 
marked  as  3.  At  this  point  the  magnetic  field  has  a  value  B  -  85  Y,  the 
solar  wind  stream  has  slowed  down  to  V  «  270  km/sec,  corresponding  to  0.4 
keV  flow  energy,  while  a  broad  ion  distribution  is  observed  with  an  equiva¬ 
lent  temperature  of  4  x  10^  °K  (i.e.,  0.4  KeV).  Notice  that  at  this  point 
the  solar  wind  has  lost  more  than  80$  of  its  flow  energy  and  its  thermal 
spread  is  comparable  to  its  flow  speed.  This  corresponds  to  the  peak  of 
the  electrostatic  wave  activity  (see  Fig.  3).  The  slowing  down  of  the 
solar  wind  continues  until  UT  12:33:20.  The  magnetic  field  exceeds  120 
Y  and  approaches  its  maximum  value  of  145  Y.  Finally,  at  UT  12:33:15  we 
note  the  beginning  of  the  magnetic  cavity,  the  suppression  of  the  wave 
activity,  and  the  appearance  of  sunward  flowing  0.5  KeV  ions.  A  summary  of 
the  key  parameters  observed  during  the  above  times  is  given  in  Table  I. 

We  note  that  there  is  a  marked  difference  in  the  nature  of  the  slowing 
down  of  the  30lar  wind  between  the  time  periods  UT  12:36:20  -  UT  12:35:27 
and  UT  12:35:27  -  UT  12:34:17.  In  the  former  case  we  point  out  that  the 
low  frequency  electrostatic  wave  activity  is  intense  and  that  there  is 
little  change  in  the  z-component  of  the  solar  wind  velocity.  In  the  latter 
time  period,  the  electrostatic  noise  has  weakened  considerably  and  there  is 
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a  substantial  increase  in  V 

pz 

detail  in  Section  IV. 


These  issues  will  be  discussed  in  more 


III.  COLLISIONLESS  MOMENTUM  COUPLING:  THEORETICAL  CONCEPTS 


Prior  to  discussing  the  interpretation  of  the  above  data  with  respect 
to  the  physics  of  momentum  coupling,  we  present  a  brief  review  of  the 
various  coupling  processes.  We  consider  the  following  momentum  equation  in 
the  x  direction  for  a  solar  wind  proton  (i.e.,  radial  direction 
perpendicular  to  the  ambient  magnetic  field  B  =■  B  e^.). 


dV 

.  Px 
dt 


eE 


a  v 

p  pz 


-  V' 


>(V 


px 


-  V 


(1 ) 


where  the  subscript  p  refers  to  solar  wind  protons,  =  eB/m^c,  v*  is  an 

anomalous  ion-ion  collision  frequency,  e is  the  streaming  barium 

velocity,  e  is  the  charge,  and  m  is  the  mass.  The  first  term  on  the  RHS  of 

(1)  arises  from  a  laminar  electric  field  usually  associated  with  the 

leading  edge  of  the  magnetic  compression;  it  acts  to  accelerate  barium  ions 

and  to  slow  down  solar  wind  protons.  The  second  term  is  the  magnetic  force 

which  is  associated  with  Larmor  coupling.  The  final  term  corresponds  to 

turbulent  "pick  up"  of  the  solar  wind  and  arises  because  of  plasma 

instabilities.  For  plasma  turbulence  such  that  v*  >  it  is  clear  that 

turbulent  coupling  can  dominate  over  Larmor  coupling.  Also,  notice  that 

the  force  Q  V  associated  with  Larmor  coupling  is  proportional  to  the 
p  pz 

value  of  V  and  will  be  very  weak  as  long  as  V  »  0.  The  time  evolution 
pz  -  pz 

of  Vpz  is  given  by 
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eE 
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n  v 
p  px 


(2) 


m  me  r 

In  deriving  (2)  we  assume  E  p=*  V  8/c  which  is  the  motional  electric  field 

z  px 

of  the  solar  wind. 

The  subject  of  the  appropriate  value  of  v*  and  the  dominant 
instability  that  drives  it  has  been  extensively  studied.  We  refer  the 
interested  reader  to  Lampe  et  al.  (1975)  and  summarize  only  the  key 
conclusions.  The  counterstreaming  between  the  barium  ions  and  solar  wind 
protons  generates  a  local  velocity  distribution  function  such  as  shown  in 
Fig.  5.  Figure  5  is  drawn  for  convenience  in  the  solar  wind  reference 
frame.  For  singly  ionized  barium  the  electrons  have  a  relative 


velocity  V 

VpX(nb/ne)  where  nh  is  the  density  of  the  barium  and  ne 


Vg  e^  with  respect  to  the  solar  wind  protons  given  by  Vg 


nK  +  n 


P* 
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The  dispersion  equation  for  this  situation  is  given  by  (Papadopoulos 
et  al.,  1971 ) 


u).  u>  -  kV  cos  6 

■><-•  »  ■  fr  Z'[-T7f - ) 

k  vw  b 
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+  — 2_  Z 

,22  '■kv  J 
k  v  p 
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♦  1  ♦-?  [l 
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c  k 


2  j /p  i /2 

where  u>  »  (4irn  e  /m  1  is  the  plasma  frequency  and  v  ■  (T  /m  )  is 
a  a  a  a  a  a 

the  thermal  speed  of  species  a  (e:  electron,  b:  barium;  p:  proton),  k  = 


k  e  +  k  e  ,  9 
XX  z  z 


(k_/k  ) ,  V 
z  x  9 


V  e  ,  a 
px  X*  e 


eB/mec  is  the  electron 


cyclotron  frequency,  0g  »  8irneTe/B-  and  Z'(^)  »  -  2(1  +  t;  Z (r))  with  Z  the 

plasma  dispersion  function.  In  writing  (3)  we  have  assumed  k  •  3  »  0 

(i.e.,  ky  =  0)  so  that  we  are  only  considering  flute  modes,  and  have 

assumed  that  the  ions  are  unmagnetized  (valid  for  u  >  £i  .  and  kp  .  >>  1 

P ,  b  p ,  b 

where  p  .  is  the  mean  ion  Larmor  radius  of  the  protons  and  barium  ions, 
p,  b 

respectively)  and  the  electrons  are  magnetized.  We  comment  that  retaining 
a  finite  ky  can  generate  the  modified  two  stream  instability  via  Vg 
(McBride  et  al.,  1972).  However,  this  instability  produces  little  momentum 
coupling  between  the  counterstreaming  ions;  it  primarily  heats  electrons 
and  generates  electron  tails  parallel  to  B. 

We  can  simplify  (3)  by  assuming  cold  ions,  i.e.,  v  <<  w/k  -  V  cos  0 

-  D  px 

and  Vp  <<  u)/k.  In  this  limit  Z (5)  =  -  1/5  -  l /2zJ  and  (3)  can  be  written 
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where  the  proton  lower  hybrid  frequency  is 
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137  then  a  <  1  and  the  small  term  in 


Notice  that  as  long  as  n./n  <  m./m  *  137  then  a  <  1  and  the  small  term  in 

b  p  0  p 

the  dispersion  relation  is  due  to  barium  ions.  This  implies  that 
when  n^/n  <  137  the  excited  waves  will  be  proton,  rather  than  barium, 
lower  hybrid  waves. 

2  2  2 

For  a  <  1  and  ioq  <<  c  k  ,  the  most  unstable  waves  have  the  following 
approximate  frequency  and  wavenumber 


“V  3  “Hp 


Y  -  0.69  a1/3uit 


k  cos  9  *  /V  .  (6c) 

Hp  px 

2  2  2 

On  the  other  hand,  for  values  of  <dq/ c  k  sufficiently  large,  the  modes  can 
be  stabilized  because  of  electromagnetic  effects  (Papadopoulos  et  al., 
1971;  Lampe  et  al.,  1975).  The  criterion  for  instability  is  given  by 

0M  „<!!&(,  .  „,/3)3/2(,  •  6J,/2  -  cos  V  (7) 

px“e 

For  the  parameters  of  interest,  we  note  that  a  <<  1  and  Bg  <<  1  and  (7)  can 
be  rewritten  as 

C03  0  <  (l  +  1  (8) 

px  p 

where  VAp  =  B/(4irnpmp)  is  the  local  Alfven  velocity  associated  with  the 
protons.  Notice  that  the  angle  0Q  separates  the  angular  region  of  unstable 
modes  from  stable  modes  !.n  the  plane  perpendicular  to  the  magnetic  field 
(Fig.  6).  For  values  of  the  RHS  of  (8)  which  are  comparable  to  or  larger 
than  unity,  the  entire  k-space  plane  is  unstable,  and  as  shown  in  Lampe  et 
al.  (1975),  complete  ion-ion.  momentum  coupling  accompanies  the  interaction. 
The  instability  weakens  substantially  when  the  stable  region  around 
increases  (i.e.,  0Q  increases)  leaving  only  weak  off  angle  modes  unstable. 

To  better  illustrate  the  linear  properties  of  the  counterstreaming 
ion-ion  instability,  we  present  numerical  solutions  of  (3)  for  parameters 
relevant  to  the  AMPTE  release  of  Dec.  27,  1984.  In  Figs.  7  and  8  we  take 
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mb/mp  =  137,  vp/Vpx  =  0.3,  vb/Vpx  =  0.01,  and  Bg  =  0.0.  The  most 
significant  variations  are  for  nb/np  and  VAp^px  so  we  present  results  for 
a  range  of  values  for  these  parameters.  In  Figs.  7  and  8  we  plot  T^/oo^p  vs 
VAp/vpx  for  nb/np  =  °'5,  2'°»  5'°  and  10,0  (or  a  ~  3*7  x  10-3,  1.5  x  IQ'2, 
3.7  x  10"2,  and  7.4  x  10-2,  respectively).  In  Fig.  7  we  consider  0  =  0°, 
while  in  Fig.  8  we  consider  6  *  60°.  Here  ym  is  the  maximum  growth  rate  as 
a  function  of  k.  We  have  not  plotted  oi^  or  k  but  note  that  (6a)  and  (6c) 
are  in  reasonable  agreement  with  the  numerical  values.  In  Fig.  7  we  note 
the  following.  First,  in  the  limit  ^Ap^px  >>  1,  the  growth  rate 
asymptotes  to  its  maximum  value  which  is  approximately  given  by  (6b)  [e.g., 
for  nb/np  =  0.5  (a  =  3.7  x  10~3)  we  obtain  TM/wHp  ■  0.1 08].  Also,  in  this 
limit  we  note  that  the  maximum  growth  rate  increases  as  nb/np  increases. 
Second,  as  VA p^px  decreases  the  growth  rate  decreases;  for  sufficiently 
small  values  of  vAp/Vpx  tJie  m°des  are  stabilized  because  of  electromagnetic 
effects.  However,  note  that  the  critical  value  of  vAp/Vpx  which  stabilizes 
the  modes  decreases  with  decreasing  nb/np  (e.g.,  for  VAp/Vpx  =  2  the  mode 
is  stable  for  nb/np  =  10.0  but  unstable  for  nb/np  =  0.5).  This  is  consis¬ 
tent  with  (8)  which,  for  0  =  0,  can  be  written  as  V„_/V„v  >  (1  +  n, /n  ). 

rtp  px  bp 

We  also  note  that  the  critical  value  of  VAp/Vpx  Predicted  by  (8)  is 
somewhat  greater  than  is  found  from  Fig.  7.  That  is,  for  hb/np  »  0.5,  2.0, 
5.0,  and  10.0,  (8)  predicts  Y  =  0  for  VAp/Vpx  »  1.5,  3-0,  6.0,  and  11.0, 
respectively;  however,  Fig.  7  shows  that  Y  =  0.01  ojHp  (*  0)  for  VAP/VPX  " 
1.2,  2.0,  3.6,  and  5.5.  The  reason  for  this  is  that  (8)  is  based  upon  the 
assumption  of  co1  protons  (u  >>  kv^)  but  for  the  parameters  used  in  Fig.  7 
(i.e.,  Vp/VpX  =  0.3)  this  assumption  breaks  down  as  V^p^px  decreases  and 
thermal  effects  allow  the  modes  to  grow  in  the  stable  "cold"  plasma  regime. 

In  Fig.  8  we  plot  TM/coHp  vs.  vAp/Vpx  for  the  3ame  Parameters  as  in 
Fig.  7  but  we  consider  0  =  60°.  The  basic  features  of  Fig.  8  are  the  same 
as  Fig.  7.  However,  two  points  are  worth  mentioning.  First,  and  most 
important,  unstable  modes  exist  for  values  of  ^Ap^px  that  lead  to  stable 
modes  in  the  case  of  0  •  0.  In  fact,  for  all  values  of  nb/np  considered, 
strong  growth  (i.e.,  Y  >  0.01  ojh  )  exists  for  VAp/vpx  >  2’  and  in  the  case 
of  nb/np  =  0.5  exists  for  vAp/Vpx  >  0>5.  Again,  this  behavior  is 
consistent  with  (8).  Second,  we  note  that  for  the  same  values  of  nb/np, 
the  maximum  growth  rates  are  smaller  for  0  =  60°  than  0  =0°.  This  appears 
to  be  inconsistent  with  (6b)  which  indicates  the  maximum  growth  rate  does 
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not  depend  upon  angle.  The  reason  for  the  discrepancy  is  the  effect  of 
thermal  protons  which  have  a  stabilizing  influence  because  of  ion  Landau 
damping.  Finally,  we  note  that  we  have  chosen  9  =  60°  because  nonlinear 
studies  have  demonstrated  that  if  the  angular  region  of  unstable  waves 
initially  satisfied  the  condition  cos  9  <  cos  9Q  *  0.5  (i.e.,  eQ  =  60°) 
then  the  entire  k-plar.e  subsequently  becomes  unstable  resulting  in  complete 
momentum  coupling  between  the  ion  streams.  If  cos  9  >  0.5  then  the 
situation,  although  mildly  unstable,  is  quickly  stabilized  by  finite 
temperature  effects  with  little  momentum  coupling. 
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IV.  COMPARISON  BETWEEN  THEORY  AND  OBSERVATIONS 

We  now  proceed  to  analyze  the  observations  with  the  theoretical  models 
described  in  Section  III.  Our  approach  is  to  examine  in  detail  the  key 
times  marked  as  1-4  in  Figs.  2  and  3.  Table  1  gives  the  values  of  the 
important  plasma  parameters.  We  note  that  all  of  the  data  are  not  well 
known  and  we  have  made  estimates  of  some  values.  Also,  in  comparing  the 
theoretical  conditions  for  instability  with  experimentally  observed 
parameters,  one  does  not  expect  to  find  "grossly"  unstable  conditions  since 
the  turbulence  observed  is  generally  in  the  nonlinear  regime  (i.e.,  near 
marginal  stability). 

The  slowing  down  of  the  solar  wind  begins  at  the  point  marked  1  in 
Figs.  2  and  3-  From  Table  1  we  see  that  n^/np  -  0.2  and  VAp/Up  =  °-4'  from 
Fig.  8  note  that  these  conditions  are  marginally  stable  for  the  ion-ion 
instability.  Therefore  the  beginning  of  the  interaction  is  consistent  with 
the  minimum  condition  for  a  momentum  coupling  cross  field,  counterstreaming 
proton-barium  instability.  The  theoretically  expected  electrostatic  waves 
cover  the  range  between  f„  «  10  Hz  and  foa  »  560  Hz,  with  most  of  the 

tip 

energy  confined  in  the  15-40  Hz  region.  At  the  time  marked  2  (UT  12:35:15) 
the  magnetic  field  compression  starts  and  there  is  an  attendant  increase  in 
VAp/Vpx.  Thus,  this  leads  to  conditions  more  favorable  to  instability  and 
we  expect  the  entire  k-plane  to  become  unstable  leading  to  strong  momentum 
coupling.  Figures  2  and  3  seem  to  confirm  this.  The  range  of  the  unstable 
waves  corresponds  to  22  Hz  -  1.6  kHz  with  the  most  of  the  energy  in  the  22 
-  60  Hz  region.  Notice  that  between  times  marked  1  -  3.  the  value  of  Vn_ 

pz. 

is  relatively  unchanged  and  remains  close  to  zero,  while  Vnv  is  reduced 
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sharply.  The  ion-ion  driven  interaction  seems  to  terminate  at  the  time 

marked  3  (UT  12:34:27).  At  this  time  the  solar  wind  speed  is  |V  |  *  250 

km/sec;  the  protons  have  lost  more  than  80?  of  their  initial  energy.  We 

comment  that  the  proton  temperature  is  of  the  order  5  x  10^  °K  so  that  the 

proton  thermal  speed  is  Vp  *  200  km/sec  and  the  wave  modes  are  subsequently 

suppressed  because  of  proton  Landau  damping.  This  is  consistent  with  the 

wave  measurements  shown  in  Fig.  4.  Following  the  time  marked  3  (UT 

12:34:27)  till  4  (UT  12:33:50)  the  value  of  Vnv  continues  to  decrease  to 

almost  zero.  However,  the  data  are  indicative  of  a  different  interaction. 

The  slowing  down  is  characterized  by  very  weak  electrostatic  activity  and 

most  important,  by  an  increase  in  the  value  of  V  which  reaches  250  km/sec 

pz 

when  VpX  approaches  zero.  This  is  the  type  of  interaction  expected  from 
Larmor  coupling  described  by  (1)  and  (2)  for  v*  <  It  is  basically  a 
gyration  of  the  protons  about  the  magnetic  field  which  is  piled  up  in  the 
front  of  the  barium  cloud.  This  stage  is  followed  by  entry  into  the 
magnetic  cavity.  We  can  associate  the  observed  sunward  flux  of  0.5  keV 
protons  with  the  thermal  expansion  of  the  protons  when  V  became  small. 

px 

The  detailed  spectrum  presented  in  Fig.  4  allows  a  further  comparison 
with  theoretical  concepts.  As  noted  before,  the  instability  saturates  by 
trapping.  For  the  proton-barium  situation  and  with  a  <<  1  the  phase 
velocity  Vph  of  the  unstable  waves  lies  near  the  barium  flow  velocity  (see 
Fig.  6).  In  the  solar  wind  reference  frame  we  note  that  (Lampe  et  al., 
1975) 


V  (l  -2* 
pxv 


Waves  growing  with  this  phase  velocity  will  trap  protons  when  their 
potential  is  of  the  order  of 
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the  condition  for  proton  trapping  (10a)  will  be  reached  first  and  will 
saturate  the  instability.  Since  in  our  case  (11)  is  satisfied,  we  expect 
that  the  instability  will  saturate  by  trapping  the  protons.  This  is 
similar  to  the  saturation  of  the  Buneman  instability  studied  by  Davidson  et 
al.  (1970)  which  saturates  by  trapping  electrons  rather  than  protons, 
despite  the  fact  that  Vph  is  near  the  proton  beam.  The  amplitude  of  the 
electric  field  required  to  trap  the  protons  is 
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(12) 


For  the  parameters  given  in  Table  I  (i.e.,  Vpx  =  250  -  400  km/ sec  and 

=  107  -  145  sec-1)  we  find  that  (12)  predicts  E  =  67  -  145  mV/m. 
Gurnett  et  al.  (1985)  report  values  of  E  =*  1 40  mV/m  while  noting  that  the 
10  Hz  channel  was  saturated.  Thus,  these  field  values  are  consistent  with 
the  saturation  by  trapping.  Note  that  the  peak  spectral  density  in  Fig.  5 
is  <  20  Hz  which  is  consistent  with  theoretical  values.  The  value  of  v* 
can  be  estimated  on  the  basis  of  the  wave  spectral  energy  density  S(ui)  in 
the  lower  hybrid  region,  as  given  in  Fig.  4.  It  is  approximately  given  by 
(Papadopoulos,  1977) 
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where  D  is  the  diffusion  coefficient  in  velocity  space.  For  the  measured 
values  at  12:34:27  and  S(<a  *  u)u  )  given  by  Fig.  4  we  find  v*  »  1 2  -  20  so 

tip 

that  v*  »  Q  ,  fl. . 
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V.  SUMMAHY  AND  CONCLUSIONS 

We  have  presented  an  analysis  of  previously  reported  observations  on 
the  30lar  wind-barium  interaction  associated  with  the  AMPTE  artificial 
comet  release  of  Dec.  27,  1984.  Based  on  these  results  we  have  argued  that 
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the  solar  wind  couples  momentum  (and  energy)  to  the  barium  ions  through 
both  laminar  and  turbulent  processes  [see  (1)].  The  laminar  forces  acting 
on  the  particles  are  the  laminar  electric  and  magnetic  fields;  the 
turbulent  forces  are  associated  with  the  intense  electrostatic  wave 
activity.  This  wave  activity  has  been  shown  to  be  caused  by  a  cross-field 
proton-barium  ion  streaming  instability.  The  observed  wave  frequencies  and 
saturated  amplitudes  are  consistent  with  our  theoretical  analysis. 

The  following  picture  emerges.  After  the  barium  is  released  from  the 
canisters,  it  expands  outward  and  is  photoionized.  The  expanding  barium 
cloud  forms  a  compressed  density  and  magnetic  field  region  on  the  sunward 
side  of  the  expansion.  As  the  solar  wind  protons  stream  into  this  region 
they  first  interact  with  the  barium  ions  to  generate  relatively  weak  off- 
angle  electrostatic  turbulence.  This  occurs  when  nb/np  <  1  and 
Vap/VpX  <  1 .  This  turbulence  acts  to  couple  the  solar  wind  protons  and 
barium  ions,  and  the  protons  slow  down  (see  in  Fig.  1  between  marks  1  and 
2).  As  the  protons  move  deeper  into  barium  ion  shells  they  encounter  the 

compressed  magnetic  field  region  (which  increases  the  local  proton 

* 

Alfven  speed)  and  a  more  dense  barium  ion  region  (which  increases  nb/np). 
These  two  factors  allow  stronger  wave  growth  to  occur;  this  corresponds  to 
the  intense  wave  activity  observed  between  the  times  marked  2  and  3  on  Fig. 
3.  This  turbulence  causes  the  solar  wind  protons  to  slow  down  even  more, 
and  also  produces  proton  and  barium  ion  heating.  Finally,  the  protons 
undergo  a  gyration  about  the  compressed  magnetic  field  which  is  indicated 
by  a  decrease  in  Vpx  and  an  increase  in  Vpz.  We  expect  a  similar  type  of 
interaction  to  occur  in  the  solar  wind-comet  interaction.  This  topic  is 
currently  under  study  and  will  be  reported  elsewhere. 

There  are  two  more  points  that  we  would  like  to  address.  The  first  is 
about  the  electron  heating  and  electron  acceleration.  The  observed 
electron  heating  seems  consistent  with  adiabatic  heating.  However,  in 
addition  to  the  local  heating,  the  UK  spacecraft  observed  electron  fluxes 
with  energy  larger  than  100  eV,  i.e.,  the  period  between  our  marks  1-3. 
This  is  the  period  during  which  the  lower  hybrid  instability  was 
operating.  For  the  flute  mode  (ky  «  0)  instability  discussed  above,  the 
electrons  are  adiabatic.  However,  field  aligned  suprathermal  electron 
tails  can  be  produced  by  the  non-flute  modes  (k  *  0),  corresponding  to  the 


electron-ion  modes  mentioned  in  Section  III.  These  modes  have  frequencies 
typically  5~6  fHp  and  saturate  at  a  lower  level.  The  details  of  this 
interaction  will  be  discussed  elsewhere.  It  is  sufficient  here  to  note 
that  the  existence  of  wave  frequencies  in  the  200  Hz  to  1  kHz  range  is 
consistent  with  the  model  and  with  the  electron  fluxes  observed  by  the  UK 
spacecraft.  Second,  although  the  solar  wind  barium  interaction  results  in 
complete  momentum  coupling,  the  barium  has  not  been  picked  up  by  the  solar 
wind  during  the  3-J*  min  of  the  measurements.  This  i3  due  to  the  large 
barium  mass  and  density  which  would  require  times  of  the  order  of  8-10  min 
to  be  picked  up.  This  should  be  contrasted  with  the  lithium  releases  for 
which  pick  up  occured  at  much  shorter  times. 

The  AMPTE  data  presented  supports  both  the  laminar  and  turbulent 
coupling  mechanisms  that  have  been  proposed  for  debris-air  coupling 
following  a  HANE  (Longmire,  1963;  Lampe  et  al.,  1975;  Goodrich  et  al., 
1985).  To  our  knowledge,  this  is  the  only  experimental  data  which 
unambiguously  demonstrates  the  various  coupling  mechanisms,  and  lends 
strong  support  to  the  HANE  coupling  models  developed  at  NRL.  Furthermore, 
it  is  indicative  of  the  type  of  data  that  would  be  valuable  in 
understanding  the  DNA/NRU  early  time  experiment;  well  resolved  data  on  the 
evolution  of  the  magnetic  field,  particles,  and  wave  activity.  Finally, 
the  good  agreement  between  theory  and  data  suggests  that  other  data 
involving  the  coupling  of  oounterstrearaing  plasmas  be  examined; 
specifically,  the  recent  data  obtained  from  the  space  missions  to  Halley's 
comet  should  be  investigated  with  regard  to  the  physics  of  HANEs. 
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TABLE  I 

Plasma  Parameters  for  the  27  Dec.  1 98^4  AMPTE  Release 
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Particle,  field,  and  electric  field  wave  data  (from  Gurnett  et 


al.,  1985).  The  barium  ion  density  is  based  upon  emissions  at 
the  plasma  frequency.  Note  the  intense,  low  frequency  (f  -  30 
Hz  -  1  kHz)  electrostatic  waves  during  between  the  times  marked 
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